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Galaxy clusters are the largest and most massive gravitationally bound structures known in the 
Universe. Cosmic-Ray (CR) hadrons accelerated at structure formation shocks and injected by 
galaxies, are confined in galaxy clusters where they accumulate for cosmological times. The pres¬ 
ence of diffuse synchrotron radio emission in several clusters proves the existence of high-energy 
electrons, and magnetic fields. However, a direct proof of CR proton acceleration is missing. The 
presence of CR protons can be probe through the diffuse 7 -ray emission induced by their hadronic 
interaction with the Intra-Cluster Medium (ICM). The Perseus cluster, a nearby cool-core cluster, 
has been identified to be among the best candidates to detect such emission. We present here the 
results of a very deep observation of the Perseus cluster with the MAGIC telescopes, accumu¬ 
lating about 250 hours of data from 2009 to 2014. No evidence of large-scale very-high-energy 
7 -ray emission from CR-ICM interactions has been detected. The derived flux upper limits in the 
TeV regime allow us to put stringent constraints on the physics of cluster CRs, in particular on 
the CR-to-thermal pressure, the CR acceleration efficiency at formation shocks and the magnetic 
held of the central cluster region. 
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1. Introduction 


Clusters of galaxies represent the latest stage of the structure formation, produced by merging 
small group of galaxies and gas accretion. Their study is a powerful cosmological tools to test 
the evolution model of the Universe ||I|]. Cosmic-ray (CR) protons can be accelerated by structure 
formation shocks, and outflows from galaxies and active galactic nuclei (AGNs) of the cluster, and 
accumulate in it for cosmological times. These CR protons can interact hadronically with the pro¬ 
tons of the intra-cluster medium (ICM), a hot thermal plasma with ksT ~keV, and generate pions. 
While TT^’s decay to secondary electrons and neutrinos, the Tl'^’s decay directly to high-energy 7 
rays. Despite many observational efforts, diffuse 7 -ray emission from clusters has remained unde¬ 
tected. 

Non-thermal emission is observed at radio frequencies in many galaxy clusters in the form 
of diffuse synchrotron radiation [^. This probes the presence of relativistic CR electrons and 
magnetic fields in the cluster environment. But, a proof for CR-proton acceleration has yet to be 
found. The observed diffuse radio emission in clusters can be divided in three main categories: 
peripheral radio relics, central radio halos and mini-halos |Q ||]. While radio relics are though to 
be related to merger shocks, the origin of radio halos and mini-halos have been historically debated 
between re-acceleration and hadronic models. In the re-acceleration model, a seed population of 
CR electrons is re-accelerated by turbulence, while in the hadronic scenario the radio-emitting 
electrons are secondaries produced by CR protons interacting with the ICM. Actually, the presence 
and role of CR protons in clusters can be probed directly only thr'ough the 7 -rays and neutrinos 
induced by hadronic interactions. The high-energy astronomy window is then crucial to understand 
non-thermal phenomena in galaxy clusters. 

The Perseus cluster of galaxies (Abell 426) is a relaxed, cool-core cluster located at a distance 
of about 78Mpc (z=0.018). It hosts the brightest thermal X-ray emission from the ICM and 
a very luminous radio mini-halo ii- The high ICM density at the centre implies a high target 
density for CR hadronic interactions. Therefore, Perseus is the best cluster where to search for CR- 
induced 7 -ray emission (see SB1I3I [TT| ] for a detailed argumentation). This cluster also hosts two 
7 -ray bright AGNs: NGC 1275, the central dominant galaxy of the cluster, and IC 310, a peculiar 
object that could be an intermediate state between a BL Lac and a radio galaxy. Both were detected 
with Fermi-LAT |I|] and MAGIC [0, IH]. 

Since 2008, the Perseus cluster is intensively observed by the MAGIC telescopes. Here, we 
present the results obtained with ~250 hours of effective observation time taken in stereoscopic 
mode from 2009 to 2014, where constraints on the CR population in the cluster are derived. 


2. MAGIC observations and results 

MAGIC is a system of two 17 m diameter Imaging Atmospheric Cherenkov Telescopes lo¬ 
cated on the Canary island of La Palma. It observes the 7 -ray sky from ~50 GeV to more than 
50 TeV. The observation of the Perseus cluster started in 2008 with ~24 h of observation with a 
single telescope, that did not allow any detection 0. Since 2009, MAGIC operates in stereoscopic 
mode, providing a much better sensitivity [p^. Here, we combine all the Perseus-cluster stereo 
observations carried out from October 2009 to November 2014. The observations of a first period. 
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Figure 1: Perseus cluster sky maps for an energy threshold of 250 GeV (left-hand panels) and 1 TeV (right- 
hand panels) obtained from ~250 h of MAGIC observation. Top panels show the relative flux (colour code, 
expressed in signal-to-background ratio) and the excess significance (contour lines starting from 4(7 with 
steps of 2a). Lower panels show the significance maps where the signal from IC310 is included in the 
background model. Symbols indicate the positions of the three brightest radio galaxies of the cluster. 


2009-2011, taken solely during dark time at low zenith angle (from 12° to 36°), led to the discov¬ 
ery of IC 310 [ 1I4| ] and NGC 1275 115]. The observations of the period 2012-2014 were carried out 
after an upgrade of the telescopes [ l^j. They were taken with relaxed criteria including Moonlight- 
night and a larger zenith-angle coverage from 12° to 60°. After data-quality selection, the final 
sample in the two periods considered amounted to an effective time of 85 h and 168 h, respectively. 


The standard MAGIC Analysis and Reconstruction Software (MARS) is used to analyze the 
data. For the 2009-2011 data sample, we use the same analysis than in our previous study on 
cluster CR-physics ||^. For the 2012-2014 data sample, to properly handle the data taken under 
Moonlight we use 33% higher image cleaning levels than the standard levels reported in [jT^. This 
increases the energy threshold but does not affect the performance at high energy. 


The 253 h of stereo observation, combined at the latest stage of the analysis, provide the deep¬ 
est view of the Perseus cluster at very high energy. The right- and left-upper panels of Figure show 
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the relative-flux (i.e., signal-to-background ratio) sky maps for an energy thi'eshold of 250 GeV and 
1 TeV, respectively. A clear signal is detected from the two previously discovered AGNs. The bright 
and hal'd source IC 310 is visible in both maps with high significance and could mask smaller sig¬ 
nals. In order to search for weak emissions, we included the point-like emission from IC 310 in our 
background model. The lower panels of Figure |T] show the resulting significance sky maps. Above 
250 GeV, NGC 1275 is detected with a significance >8a. The shape of the detected signal is in 
perfect agreement with a point-like object such as expected for an AGN. Above 1 TeV, however, no 
source other than IC 310 is detected. For both energies, there is no sign of diffuse y-ray structures 
inside the cluster. 

The average energy spectrum of NGC 1275 obtained with the full 2009-2014 data set is shown 
in Figure ^ together with the previously reported spectra from the first two years of observation 
[11^. The spectrum of this analysis starts at higher energy because the data include Moonlight 
and large-zenith-angle observations. The statistical precision is largely improved, extending the 
spectrum coverage with a data point around 880 GeV. The spectmm between 90 GeV and 1200 GeV 
can be described by a simple power law 

E 

200 GeV 

with a photon index F = 3.6 ±0.2^,ar 3=0.2^^,^, and a normalization constant at 200 GeV of /o = 

(2.1±0.2,f«f ±0.3,y,,) X 10^“cm-2s-iTeV-i. 




3. Search for cosmic-ray induced emission 


The y-ray emission from the central galaxy NGC 1275 is consistent with a point-like source 
and no diffuse component is observed. The measured flux is much larger than what is expected 
from the CR-induced emission, and it could be outshining the latter. However, the NGC 1275 AGN 
specttum is very steep and no signal is detected above about 1 TeV, while models expect the CR- 
induced spectrum to be a flat power-law with no cutoff in the MAGIC energy range |||]. Therefore, 
we use the energies above 1 TeV to search for the possible diffuse CR-induced component. 

As the spatial shape of the CR-induced emission is not known, three different models are used 
as template: i) the isobaric model assuming a constant CR-to-thermal pressure Acr = FcR/Rth [p0|], 
ii) the semi-analytical model for CRs derived from hydrodynamical simulations of clusters [S], and 


hi) the extended hadronic model of [21], where the possibility of CR propagation out of the cluster 


core is considered, resulting in a significantly flatter CR profile. The y-ray surface brightness of 
these three models is shown in Figure |3[ The fraction of signal expected within a circular region 
of a given radius 6 from the cluster centre and the MAGIC Point Spread Function (PSF) above 
630 GeV are also shown. Since the predicted CR-induced signal extension is significantly larger 
than the PSF, optimal selection cuts 6cut, different than for a point-like source, are used to detect 
the emission for each of the models, which are shown in Figure 

To derive the strongest possible constraints, we compute upper limits on the integral flux above 
a given energy that accumulates results of several energy bins. The obtained limits depend on the 
assumed spectral shape. We adopt the universal CR-proton spectrum found with hydrodynamical 
simulations [1^, that suggests a power-law momentum spectrum with a 2.2 at the energies 
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Figure 2: Spectral energy distribution of NGC 1275 averaged over different periods. The red symbols show 
our new result obtained with ^250 h of data. The red arrows are the differential flux upper limits for a 
power-law spectrum with a photon index r=3.5 (thick solid lines) and r=2.3 (thin and dashed lines). 




Figure 3: Left: Surface brightness profiles of the three tested spatial templates for the CR-induced emission 
in Perseus. Right; Cumulative fraction of signal within a given radius for different models; point-like in 
black (long-dashed line), isobaric in green, semi-analytical in red, and extended in blue. The distribution 
of the point-like model corresponds the the MAGIC PSF above 630 GeV. The thin and thick coloured solid 
lines represent the real and MAGIC-reconstructed, i.e., smeared by the PSF, signal fractions, respectively. 
The dotted lines show the fraction of these signals contained in the background region. The vertical dashed 
lines correspond to the used optimal dcut values. 


of interest here. At very high energy, the y-ray spectrum induced by pion decays should have ap¬ 
proximately the same spectral index. We convert this intrinsic spectrum into the spectrum observed 
on Earth, taking into account the y-ray absorption by the Extragalactic Background Eight (EBE) 
[22] from the Perseus cluster at z=0.018. Between 300 GeV and lOTeV, the effect of the absorption 
can be approximated by an increase (softening) of the power-law index of about 0.13 and a reduc- 
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Table 1: Integral flux upper limits within a radius of 0.15° and 1.4° from the Perseus cluster centre, between 
Efh and 10 TeV for a power-law spectrum with an index F = 2.33 [in unit of 10^^^ cm^^ s^^]. 


Eth 

model 

ON 

OFF 

OUMa 

p0.15“ 

^UL 

pi. 4“ 
^UL 


point-like 

332 

304.1 

1.4 

12.2 

12.2 

630 GeV 

isobaric 

1327 

1256 

1.4 

24.8 

65.5 


semi-analytic 

749 

681 

1.8 

26.8 

49.0 


extended 

1327 

1256 

1.4 

25.6 

124. 


point-like 

159 

157.5 

0.1 

3.84 

3.84 

1.0 TeV 

isobaric 

675 

652 

0.6 

10.7 

28.3 


semi-analytic 

369 

352 

0.6 

9.77 

17.9 


extended 

675 

652 

0.6 

11.1 

54.0 


point-like 

77 

75.9 

0.1 

2.34 

2.34 

1.6 TeV 

isobaric 

321 

317 

0.2 

5.09 

13.5 


semi-analytic 

169 

167 

0.1 

4.61 

8.43 


extended 

321 

317 

0.2 

5.26 

25.6 


Notes. ON and OFF are the number of events in the signal (ON) and background (OFF) regions, 
respectively. OuMa is the significance of the excess expressed in standard deviations. 


tion of the differential flux at 1 TeV of 17%. Above 10 TeV the absorption increases dramatically. 
A source with an intrinsic spectral index a = 2.2 would appear above 300 GeV as a power-law with 
an index F = 2.33 and a cutoff above 10 TeV. About 20% of the flux above 1 TeV and 60% above 
10 TeV is absorbed during the travel to Earth. This effect was neglected in all previous papers, 
which, therefore, overestimate the constrains induced by the flux upper limit. 

Table |T] presents the 95%-confidence-level upper limits (method from [19] with a total sys¬ 
tematic uncertainty of 30%) of the integral flux between several energy thr'esholds Efh and 10 TeV 
for different models assuming a power-law spectrum with F = 2.33. The upper limits are con¬ 
verted to the corresponding flux contained within the reference radius 0.15° and the cluster virial 
radius' /?200 — 1-4° The conversion factors depend on the surface brightness distributions, 
considering also the signal contamination in the OFF region. 


4. Interpretation and discussion 

One of the main uncertainties in modeling CR physics in clusters of galaxies is the CR spatial 
distribution, which is assessed by considering different spatial templates (isobaric, semi-analytical 
and extended). Another major uncertainty is the CR-acceleration efficiency, i.e., the fraction of 
the energy dissipated at structure formation shocks which causes particle acceleration. Our semi- 
analytical model is based on predictions by [||] derived from hydrodynamical simulations of clus¬ 
ters, which assume a maximum CR-proton acceleration efficiency (^p,max = 50%. Here, we tie this 

* The cluster virial radius is here defined with respect to a density that is 200 times the critical density of the Universe. 
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to the parameter Ay, just a multiplier for the whole formalism, that is Ay = 1 for (^p,max = 50%. 
Smaller values of Ay correspond to smaller values of Cp,max- 

The left-hand panel of Figure shows the integral 7 -ray fluxes predicted by 0 within 0.15° 
from the centre, Ay = 1, and the highest fluxes allowed by our observational upper limits. Compared 
to 0 ]> we accumulate three times more data and derive upper limits that are significantly lower. 
Ay is constrained to be < 0.56 when neglecting the EBL absorption and < 0.75 with the EBL[^^. 
This corresponds to a maximum CR-proton acceleration efficiency of 28% and 37%, respectively. 
The latter should be considered the reference value while the one without EBE absorption is given 
to allow comparison to previous results. 




0.01 0.10 1.00 

•^^^200 


Figure 4: Left. Integral flux upper limits and spectra within 0.15° for the semi-analytical model. We show 
the spectrum predicted by and the ones adapted to match our upper limits. Shown in light gray are our 
previous flux upper limits from [|^. Right. Volume-averaged Acr within a given radius R, as constrained by 
the upper limits presented in this work for the semi-analytical, extended, and isobaric models with a = 2.2. 

The constraints on the volume-averaged CR-to-thermal pressure profiles for the three CR- 
distribution models are shown in Eigure ^left. < Acr > results to be below 1-2% within 0.15° 
(« 0.11 X R 20 o)- When considering the full galaxy cluster volume up to R 200 > < ^CR > is below 
2% for the isobaric and semi-analytical models, but significantly less constrained, below 19% 
for the extended model. This latter weak constraint was expected because the CR distribution is 
significantly flatter in this case. 

The Perseus cluster hosts a bright radio mini-halo. Assuming the hadronic model of radio 
halos, i.e., synchrotron emission generated by secondary electrons produced in hadronic interac¬ 
tions between the CR and ICM, the pion-decay /-ray emission is directly linked to radio signal. 
The intensity of the synchrotron emission depends on the amount of electrons, proportional to the 
hadronically-induced / rays, and the local magnetic field. This emission can be reduced due to 
the electron energy lost through other channels, mainly inverse Compton. In the regime of very 
strong magnetic fields, the electrons lose all their energy through synchrotron and we can derive 
the expected /-ray emission. Eor assumptions of CR spectral index a < 2.1, the expected /-ray 
emission is above our upper limits. Eor softer CR spectrum, our /-ray upper limits can be turned 
into lower limits on the magnetic field. This constraint, depending on the photon field considered 
for the inverse Compton scattering, is still under study. 
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